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• Procede de determination des parametres optimaux d'une acquisition de: 

radiographie . 

La presente invention a pour objet un procede de determination des . . 
5 parametres optimaux d'une acquisition de radiographie. Le procede de - 
i'invention est de preference automatique. L'inverition a egalement pour 
objet un procede de reglage d'une installation de radiologie, notamment un 
procede de reglage de la haute tension d'alimentation d'un tube a rayons X 
de cette installation, ainsi que du courant de ce tube, afin de realiser cette 

10 acquisition avec cette installation. Le but de Tirivention est de maTtriser une , 
dose de rayonnement 6m\s, et d'augmenter le contraste d'une image 
acquise avec une telle installation, afin qu'elle revele au mieux des structures 
qu'on cherche a examiner. 

Le procede de {'invention, est destine a etre mis en oeuvre avant 

15 chaque acquisition, en temps reel, a partir d'une image test. ^ - v 

Une installation a rayons X requiert que des parametres d'exposition .;• 
corrects, essentiellement la haute tension kVp, I'integrale du courant debite 
mAs. et le pouvoir de filtrage de filtres interposes, soient estimes " - 
correctement, de maniere a assurer une penetration optimale du ; 

20 rayonnement dans les tissus etudies et une bonne qualite d'image. Cesv; 
parametres dependent physiquement de Tepaisseur radiologique des zones i;:- 
imagees. Un corps d'un patient presente une distribution d'epaisseurs;\.r<^ 
exprimees en termes d'epaisseurs equivalentes, en regard de chacun des 
elements d'image, pixels,, d'un detecteur de Tinstallation. Le corps peut par 

25 ailleurs etre represents par une epaisseur moyenne, dite EPT moyen, EFT 
" mean pour Equivaleht Patient Thickness - epaisseur equivalente du patient: 
Cette distribution est encore mieux representee par la dynamique patient, 
notee AEPT. qui correspond a la variation d'epaisseur equivalente des tissus 
d'interet. 

30 Jusqu'a present, dans les installations, seul PERT moyen etait pris en 

consideration. Dans les debuts de la radiologie, Tepaisseur equivalente 
moyenne du patient, I'EPT moyen, etait deduite du poids de ce patient. Par 
la suite, cette epaisseur moyenne a ete mesuree d'une fagon plus realiste, 
en prelevant une image test du patient, avec des conditions donnees 

35 arbitraires de fonctionnement de Tinstallation de radiologie, et en mesurant 
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sur rimage test * airisi relevee' des caracteristiques* d'elements' d'image 
representatifs de cette epaisseur moyenne. Pour fixer les idees, et bien que 
ceci ne soit pas une limitation, des epaisseurs moyennes pouvaient ainsi 
etre mesurees entre 4 ou 5 centimetres et une quarantaine de centimetres. 
5 En ce qui conceme la dynamique AEPT du patient, elle n'etait pas mesuree. 
Elle etait imposee arbitrairemenfa une valeur pouvant varier entre quelques 
centimetres et une vingtaine de centimetres. Cette dynamique patient etait 
' laissee a la disposition du praticieri. Ce praticien s'en servait pour regler le 
contraste des images revelees de la maniere qui lui convenait le mieux. 

10 Dans ces conditions, le reglage etait subjectif mais pas objectif. 

' La dynamique en epaisseur du corps d'un patient est la variation 
constatee/ dans les zones d'interet du patient, entre les epaisseurs 
equivalentes les plus faibles et les epaisseurs les plus fortes. Par exemple, 
pour une dynamique en -epaisseur, partant d'une epaisseur minimale (de 

15 Tordre de 3 cm) on peut alter jusqu'a une epaisseur maximum plus ou moins 
grande, dependant des tissus a imager dans le patient. Par exemple une 
dynamique utile de 5 cm en epaisseur est rencontree dans la region de 
^ I'abdomen d'un patient- (ou les tissus sont peu differencies en termes de 
densite radiologique) alors que dans d'autres regions du corps du patient 

20 une dynamique plus grande par exemple de 14 cm est rencontree. La 
dynamique de mesure doit etre adaptee respectivement de 3 a 3+5=8 cm 
dans le premier cas, ou de 3 a 3+14=17 cm dans le deuxieme. Si une 
dynamique pardefaut est retenue, typiquement une dynamique de 14 cm en 
epaisseur, et si on mesure avec un tel reglage de rinstallation une zone 

25 d'abdomen d'un patient, alors on perd beaucoup en contraste d'iniage, en 
sensibilite de detection et de revelation, puisque le signal utile pour 
I'abdomen ne sera revele qu'avec une dynamique de 5/14eme de ce qui 
serait possible. 

Linconvenient presente par ces methodes est que le contraste des 
30 images montrees est done toujours sous-optimal. En effet, en particulier 
daris le cadre des radiographies du pelvis, aucun mecanisme n'est mis en 
' oeuvre pour rejeter des zones de saturation, et aussi des zones dites de non- 
anatomie. Ceci conduit a des estimations erronees de I'epaisseur 
equivalente moyenne, EPT moyen. En outre I'absence de dispositif 
35 d'elimination des radiations de saturation force I'utilisateur a utiliser des 
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filtres de contour, et a retenir des mises en place strictes de patients. Ceci 
n'est particulierement pas facile lorsque I'image doit reveler des 
^ . angiographies a proximite du bas des poumons. 

Enfin la dynamique ne peut pas etre estimee, ce qui conduit a une 
5 gestion sous-optimale des doses et a un defaut de qualite. Dans d'autres 
cas, un defaut de reglage de la dynamique peut conduire a des saturations 
, aberrantes dans certaines images : une partie interessante dej/image se 

trouve dans une zone saturee. Dans, d'autres cas, le mauvais contraste ou 

. I- • . 

des images trop sombres peuvent conduire a un rapport -contraste a bruit 

10 inferieur de 30 % par. rapport a un contraste optima!. 

L*invention a pour objet de remedier a ces inconvenients en aglssaht 
, ^ de deux manieres. Dans une premiere maniere, le calcul ,de I'epaisseur 
moyenne qui est une information pivot pour le reglage de Tinstallation est 
; „effectue, comma dans I'etat de la technique, a partir d'une image test 

15 acquise avec une installation de radiologic fonctionnant dans des conditions 
connues de fqnctionnement. Mais dans cette image test, selon rinvention • on . - 
ecarte des elements d'image qui ne sont pas representatifs de parties 

. V. anatomiques significatives du patient. Parmi les parties non significatives on . 
.retiendra d'une part les parties saturees de Timage. Typiquement, les parties 

20 saturees correspondent a des. parties de Timage situees au-dela du bord'du - 
corps du patient. Mais de preference dans invention, on verra . qu^-|ces 
parties non significatives peuvent. correspondre a des epaisseurs de corps 
de patients inferieures a un seuil, qu*on determinera, en dega duquel on sait % 
qu'aucune structure anatomique n'a d'interet. Dans une deuxieme maniere, 

25 ayant par ailleurs determine ce seuil, on retient comme element de fixation 
de la dynamique Lin ecart "entre' ce" seuil; exprinie en" terme d'epaisseur " 
.equivalente, et Tepaisseur moyenne du corps du patient EPT moyen. On . 
dispose alors en agissant ainsi d'une grandeur representative d'une mesure 
objective de la dynamique d'affichage a retenir. 

30 Par ailleurs on montrera que ces epaisseurs de seuil peuvent elles- 

rnemes etre corrigees pour tenir compte, ou non, de certains phenomenes 
perturbateurs. de fagon a augmenter la robustesse de la determination de 
I'epaisseur equivalente moyenne ou de la dynamique. 

Uinvention a done pour objet un precede de determination des 

35 parametres optimaux d'une acquisition de radiographie, dans lequel, 
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' • - on acquiert uhe premiere image test du patient dans des conditions 
de reglage connues d'une installation de radiologie utilisee a cet effet, 

- on mesure pour ces conditions de reglage connues une epaisseur 
moyenne du patient a partir de cette premiere image test, et ^ ■ 

5 ' - on determine les parametres optimaux de I'acquisition a partir de 

cette epaisseur moyenne, ■ ' • ^ 

'" ' " caracterise en ce que . - 

' - pour mesijrer I'epaisseur moyenne a partir de cette premiere image 
test, on ecarte dans cette premiere image test des elements d'images qui ne 
10 sorit pas representatifs de parties anatomiques significatives du patient. 

- L'invention a egalement pour objet le reglage d'une^ installation de 
' radiblogie, notamment de la haiite tension appliquee entre une anode et une 
cathode d'un tube a rayon X de cette installation, ce reglage etant fonction 
de {^epaisseur moyenne d*un patient examine et etant realise de preference 
15 en temps reel, en quelques millisecondes. 

L'invention sera mieux cpmprise a !a lecture de la description qui suit 
* et a I'examen des figures qui I'accompagnent, Celles-ci ne sont presentees 
I ' "■ qu'a titre indicatif et nullement limitatif de' Tinvehtion, Les figures riionfrent : 

- Figure 1 : urie representation ischematique composite- d'un 
20 phenomene d'irradiation et de sa prise en compte, dans I'etat de la 

technique, et dans invention respectivement ; 

- Figure 2 : une representation schematique d'une partie du corps d'un 
patient dans laquelle differentes zones sont representees : des zones 
d'interet et des zones de peu d'interet ; 

25 - Figure 3: une suite d'etapes mises en oeuvre dans un procede 

prefere de reglage d'une installation radiologique selon Tinvention, 
notamment avec tous les perfectionnements dont certains sont seulement 
optionnels ; 

- Figure 4 : la representation d'un histogramme d'elements de Timage 
.30 relevee permettant de mieux comprendre le fonctionnement de Tinvention. ' 

La figure 1 montre, notamment dans sa partie superieure, une 
installation de radiologie schematique utilisable poiir mettre en oeuvre le 
procede de Tinvention. Cette installation comporte un tube T a rayons X 
emettant des rayons X RX en direction du corps C d'un patient.* Le- corps C 
35 est montre ici comme possedant un profil triangulaire. Cette representation 
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est tout a .fait artificielle mais elle permettra d'expliquer plus simplement 
rinvention, Bien entendu le corps du patient a plutot une forme ovale voire 
rectangulaire dans une coupe examinee. Le rayonnement X emis par le tube 
T est traditionnellement filtre par. des filtres formes par des lames de cuivre 
FCu et par des lames d'aluminium FAI. Ces filtres permettent de s'assurer 
que la densite spectrale des rayons X est confinee dans une bande 
passante relativement etroite. Le pouvoir de filtrage de ces filtres intervient 

.bien entendu dans le reglage de Tappareil, differents filtres pouvant etre mis 
en place a la demande. 

En dessous d'un plateau support patient non represente est situee 
une grille anti-diffusante GA superposee a un detecteur D. La grille GA 

' comporte en pratique un certain nombre de parois septales perniettant de 
s'assurer que le rayonnement qui la traverse est uniquement (en theorie) un 
rayonnement X provenant directement du tube T. La grille GA realisant 
cependant une absorption, son epaisseur est reduite, ce qui conduit a 
reduire Tefficacite de captage de diffusions Compton, formant; un 
rayonnement diffuse, et qu'on cherche a eliminer. Le detecteur D comporte 
par ailleurs un ensemble d'elements detecteurs d'un signal d'image 
corresRondant a des elements d'image, ou pixels P. Le detecteur D est en 
pratique un detecteur electronique. Dans la suite de Texpose, on assirojlera 
aux pixels le signal deiivre par les elements detecteurs situes a leur endroit. 
On pourrait toutefois aussi envisager de numeriser une image acquise sururi 
film. 

La figure 1 montre le corps C, selon une direction X et une direction 
en epaisseur, en hauteur. Mais il faut admettre que le detecteur D ainsi que 
la grille anti-diffusante GA sont evidemment a deux dimensions. Le detecteur 
D ainsi que la grille GA definissent un champ de vue FOV, pour Field of View 
en anglais. Ce champ de vue FOV s'etend sur des zones de saturation Zsat 
ainsi que sur des zones de corps, en regard du corps C qu'on separera, 
selon rinvention. en zones anatomiques Za et zones de non anatomie Zna. 

Le premier diagramme de la figure 1 represente en correspondance, 
en dessous de cette installation schematique, une epaisseur du corps C du 
patient croissant en fonction de I'abscisse X. Cette epaisseur part de zero a 
la frontiere entre la zone de saturation Zsat de gauche et la zone de non 
anatomie Zna et augmente jusqu'a I'epaisseur maximale. a I'extremite a 
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"droite du corps C. Ce diagramme en epaisseur est a la fois id un diagramme 
^ ' en epaisseur vraie et un diagramme en epaisseur equivalente, Tun se 
d^duirait de J'autre par une simple homothetie. 
• ■ ^ Les epaisseurs equivalentes sent des epaisseurs de^ tissus humains 
5 \ donnes par leurs equivalents en epaisseur de materiaux plastiques de 
qualite connue par des normes. 

Le diagramme situe en' dessous du diagramme de Tepaisseur montre 
pour cette epaisseur le signal re?u par le detecteur D. Dans les zones de 
^ saturation Zsat; le detecteur est effectivement nomnalement sature puisque 
10 ' aucun tissu n'a ete interpose sur le cheminement du rayonnement X. Le 
' • detecteur mesure une energie regue qui est depehdante de la force de 
penetration des rayons (de la durete des rayons X) et de la duree de la pose 
* ou plus simplement des milliamperes multiplies par les secondes de cette 
pose. De maniere attendue, le signal re?u decroTt en correspondance de 
15 {'epaisseur du corps C que les rayons X ont a traverser. On moritre ici 
artificiellement que ia decroissance est lineaire avec I'epaisseur. Ceci n'est 
pas vrai en theorie et en pratique du fait d'une- absorption de type 
" exponentielle. Mais cette representation simpliste permet de mieux expliquer 
rinvention. ' 
20 Si rinstallation pour les conditions connues de reglages n'est pas trop 

mal reglee pour I'epaisseur maximaie du corps C la quantite d'energie reQue 
pour les pixels de droite concernes par cette epaisseur maximaie ne sera 

pas nulle. Sinon il y a ebasage, et les conditions d'acquisition de Timage de 

test a partir desquelles on ferait toutes !es mesures seraient un peu 
25 faussees. Naturellement, au-dela a droite du bord du corps C, le signal re?u 
correspond egalement a un signal de saturation. 

Le detecteur D ou tout autre systeme equivalent d'imagerie, y compris 
la numerisation d'un cliche revele sur un film radiographique, possede une 
dynamique de revelation. En pratique, a Tendroit de chaque pixel, Tenergie 
30 regue est'mesuree par des compteurs d'echantiilonnages dont le nombre de 
positions de comptage est limite. Dans un exemple non limitatif, on a retenu 
des' compteurs a 14 positions de sorte que le signal delivre par ces 
compteurs ne peut etre que compris entre 0 et 2^^-1 c'est a dire entre 0 et 
16383 (ou 16384 en oubliant le -1). Un premier reglage simple de la 
35 dynamique du detecteur peut etre de j-egler un niveau de gris maximal, 
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correspondant aux parties blanches de Timage, pour les zones du corps en 
limite de la zone saturee a gauche, et un niveau de gris 0, correspondant au 
noir, pour Tepaisseur du corps C la plus grande. Entre les deux, ici 
artificiellement, on a montre que le signal evoluait en fonction de I'abscisse X 
5 d'une nnaniere lineaire en passant de 16384 a 0 d'un bord a Tautre du corps 
C. 

En definitive toutes ces explications expliquent, pour des conditions 
connues de reglage de rinstallation, Tacquisition d'une premiere image test 
Par exemple cette image peut etre celle representee sur la figure 2 nriontrant, 
10 tres schematiquement le pelvis, d'un patient. Sur cette image on .distingue, a 
droite et a gauche des jambes Jd et Jg, des zones de saturation Zsat de 
- . meme qu'entre ces jambes. Avant llnvention, la valeur moyenne EPT moyen f 
de I'epaisseur du corps du patient C etait mesuree en prenant en compte a -a^ 
Ja fois les zones representatives du ^corps C et a la fois les zones de . k 
15 saturation, Dans Tinvention, en particulier, on eliminera les zones Vde -^^ 
saturation, mais pas seulement ou pas exclusivement celles-la. --^l 

Si, en particulier, on determine qu'une zone d'interet Zh delimiteeiici 
par des tirets. correspond a une partie abdominale du corps du patient, oh a 
interet a augmenter le contraste dans cette zone d'interet. On notera par % 
20 ailleurs que le probleme du contraste est particulierement delicat a resoucjre Va- 
pour les zones abdominales, ou les tissus sont de fait peu differencies, etpis . 'V 
naturellement le contraste n'est pas tres bon. Y " 

Dans rinvention on a par ailleurs remarque que dans les zones dites 
de non-anatomie, typiquement representees par les bords B du corps C du t 
25 patient, il n'y avait pas d'information d'image a rechercher Dans Tinvention 
~ on a alors elimine; pour le calcul de I'epaisseur moyenne representative de ' ' * 
. conditions vraiment utiles d'acquisition, de preference, a la fois les zones de 
saturation et a la fois les zones de non-anatomie. 

Sur la figure 1, on a ainsi represents sur le diagramme de I'epaisseur, 
30 deuxieme diagramme, d'une maniere schematique, un seuil d'epaisseur S 
en depa duquel on considere que. les zones du corps examine sont des 
zones de non-anatomie. On verra plus loin comrnent ce seuil S est 
determine. Cependant, connaissant ce seuil S, on peut calculer, dans les 
conditions de I'acquisition de I'image test etudiee, a quel niveau de gris ce 
35 seuil S d'epaisseur correspond. On appelle GTH ce seuil de niveau de gris. 



1 er depot 



. . 8 ■• ■ ■ 

On peut aussi calculer son equivalent en dose regue. Pour calculer alors la 
valeur moyenne, EPT moyen, dans rinvention ; on realise 'un histogramme 
des valeurs de niveau de gris des pixels de rimage. En pratique ici, puisque 
d'une maniere Ires artificielle le corps C a un profil triangulaire, et en 
negligeant la mesure exponentielle, Thistogramme se presente sous la forme 
d'un nombre de pixels constant, quel que soit le niveau de gris, voir le 
diagramme en bas a droite de la figure 1. Get histogramme cofnporte, pour 
■les zones de saturation, un nombre tres important de pixels revelateurs d'un 
signal sature, niveaux de gris superieurs a un niveau de saturation, indiqu6 
scliematiquement a 16384 dans I'exemple. II montre par ailleurs une zone 
hachuree situee au-dessus du niveau de gris correspondant au seuil S 
d'epaisseur. 

L'idee de I'invention est de calculer I'epaisseur moyenne EPT moyen 
sur les bons pixels candidats seulement, c'est a dire sur ceux situes a droite 
de I'epaisseur seuii sur le haul de la figure 1 , et situes en dessous du niveau 
de gris correspondant dans I'histogramme. La population de i'histogramme 
ainsi reduit permet alors de calculer une epaisseur moyenne. 

En pratique, les niveaux de gris sont revelateurs des doses regues. 
Ces doses regues revelent, selon la formule I = lo exp(-Mx), que le niveau de 
rayonnement regu varie exponentlellement en fonction de I'epaisseur (x). 
Pour calculer I'epaisseur moyenne, il faut done convertir I'histogramme en 
niveaux de gris en un histogramme en epaisseur S equivalentes, en passant 
par le logarithme des niveaux de gris (ou des doses regues si on travaille en 
doses). Avec les epaisseurs equivalentes et avec le nombre de - leurs 
occurrences, on calcuie I'epaisseur moyenne. 

Celle-ci est representee ici d'une maniere tres schematique comme 
correspondant a une position intermediaire entre le niveau de gris du seuil 
d'epaisseur et le niveau de gris de I'epaisseur maximale. Cette epaisseur 
moyenne est ainsi calculee d'une maniere bien plus rigoureuse que dans 
I'etat de la technique (et on verra par la suite comment ce calcul peut lui- 
meme etre encore ameliore). Cette epaisseur moyenne est utilisee d'une 
. maniere connue pour regler rinstallation de radiologie. II suffit d'utiliser cette 
epaisseur moyenne dans un logiciel de trajectoire. Un tel logiciel de 
trajectoire comporte des moyens de determiner les parametres de reglage 
d'une installation, en fonction notamment d'une epaisseur moyenne EPT 
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moyen, d*ui?e dynamique AEPT, d'un nombre de vues souhaite, et de la 
temperature du tube T au moment de Texamen. Le calcul de trajectoire 
permet de regler Tinstallation au mieux des souhaits de Tutilisateur afin que 
le tube en fin d'experience n'atteigne pas des valeurs de temperature 
conduisant a sa deterioration. Ces logiciels de trajectoire sont connus et 
specifiques a chaque installation. 

Cependant dans rinvention. plutot que d'en rester la, avec cette 
epaisseur moyenne, on calcule egalement .la dynamique. En effet,- on peut 
admettre que les conditions.de reglage dictees par la seule connaissance de 
I'epaisseur moyenne vont correspondre a celles d'une decroissance des 
niveaux de gris, des epaisseurs nulles jusqu'aux epaisseurs les plus fortes, 
telle que representee sur la courbe C1 de la figure 1. A I'oppose, dans 
rinvention on va aussi r6gler la dynamique. Celle-ci est reglee de telle fagon 
que Tepaisseur moyenne, EPT moyen, corresponde a une proportion 
donnee de la dynamique en niveaux de gris du signal, ou en doses si on 
travaille en doses. Dans ce but d'une part, a titre de perfectionnement, plutot 
que de choisir un seuil d'epaisseur S expose jusquMci, on choisira un seuil dit 
Anatomique, max , correspondent a une zone maximale d'anatomie 
Jnteressante. Le niveau de grjs du seuil Anatomique max est inferieur au 
seuil d'epaisseur. On regie alors I'lnstallation pour que le detecteur delivre:un 
sjgnal maximal, 16384 niveaux de gris dans I'exemple, pour les epaisseurs 
correspondent a cette epaisseur de zone maximale d'anatomie. On fixe^de 
cette fagon un premier point M de reglage de la chalne de detection. D'autre 
part, un deuxieme point N est tel que, pour la valeur EPT moyen, EPTmean, 
le niveau de gris rendu par le detecteur soit egal a une proportion donnee de 
-lavdynamique. Dans un exem pie cette proportion est celle de 1/20^'^® de la ' 
dynamique. En pratique la correspondance pour EPTmean est alors de celle 
du niveau de gris maximum d'affichage multiplie par la proportion donnee, 
soit calee a 800 niveaux de gris dans I'exemple. 

Dans ces conditions, on definit la dynamique AEPT comme 
correspondant de preference a la difference en epaisseurs equivalentes 
entre I'epaisseur anatomique maximum et Tepaisseur moyenne. On aurait pu 
la faire correspondre a la difference entre I'epaisseur moyenne et I'epaisseur 
seuil. 

On obtient alors le resultat particulierement interessant, figure 2, que 
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dans la zone d'inter§t Zl, des tissus radiologiquement peu differencies sont 
alors rendus avec le cbntraste optimal pour leur discfiminatidn et pour 
i'utilisation par un pratjoien. 

La figure 3 represente une suite d'operations mises en ceuvre dans un 
dispositif de traitement d'informations de l'lnstaliation radiologique de ia 
figure 1, non represente, dans lequel toutes ces operations sont entreprises. 
La figure 3 montre une premiere operation 1 a Tissue de laquelle on a 
determine i'epaisseur seuil S en dega de laquelle on considerera que les 
zones du corps C sont des zones de nori-anatomie. Cette etape 1 est 
ensuite suivie par une etape 2, etudiee plus loin, au cours de laquelle pour 
les conditions d'acquisition de I'image test, on calcule la dose correspondant 
a cette epaisseur seuil S, ainsi que. ou de preference, le niveau de gris de 
cette epaisseur seuil. Au cours d'une etape 3. apres la prise d'image, on 
realise I'histogramme montre sur la figure 1 ainsi que I'histogramme reduit 
duquel on a enleve les elements d'image dont les niveaux de gris en 
pratique sont superieurs au seuil. 

La figure 4 montre un histogramme reel, en terme de nombre de 
pixels par niveau de gris. et non plus simpliste comme celui vu jusqu'ici. Elle 
montre que le niveau de gris du seuil S d'epaisseur peut se trouver bien au- 
dela d'une zone d'interet reel. La figure 4 montre ainsi une partie, 
correspondant a du bruit, situee, en termes d'epaisseur, entre I'epaisseur 
seuil S et une Epaisseur minimale 1 des tissus situes dans la zone d'interet 
Zl. On considere que les pixels concernes representent du bruit parce que le 
nombre de pixels concernes par chaque niveau de gris y est faible et y est 
sehsiblement constant. 

Plutot que de se satisfaire alors de Telimination des pixels dont le 
niveau de gris est superieur au niveau de gris du seuil S, I'idee 
complementaire de I'invention est egalement d'enlever, de la population qui 
va etre utilisee pour determiner I'epaisseur moyenne, ceux de ces pixels 
dont le niveau de gris est inferieur au' niveau de gris du seUil de I'epaisseur, 
mais qui par ailleurs sont proches de lui. Comme on veut eliminer les 
niveaux de gris pour lesquels les occurrences sont faibles, on etablit une 
somme cumulee des occurrences de ces pixels representatifs de bruit, partie 
hachuree de la figure 4. Dans ce but, on determine I'Epaisseur minimale 1 
correspondant a I'epaisseur d'anatomie maximale, par soustraction d'un 
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certain nombre predetermine d'occurrences cumulees de niveaux de gris. 
Cette soustraction est faite a partir de rhistogramme mesure entre ce seuil et 
cette anatomie maximale. Dans un exemple, on decide qu'on limite cette 
integrate a 1024, on decide done qu'on elimine ainsi 1024 occurrences de 
5 niveaux de gris. On pourrait toutefois prendre un riombre plus grand, par 
exenfiple 10 000. Lidee est. sur la figure 4, de rechercher une frontiere qui se 
, trouverait a gauche de la frontiere representee par le seuil S d'epaisseur 
precedemment determine. . . ^ 

. La figure 4 montre deux types de tissus. Un premier tissu T;l pour 

10 lequel .le bruit est peu important, et un deuxieme tissu T2 montre ayec des 
tirets pour lequel le bruit est tres important. On voit qu'en decidant d'enlever 
un certain nombre d'occurrences de niveaux de gris de la population de 
pixels, qui vont concourir a I'etablissement de la valeur d'epaisseur moyenne 
EPT moyen, il est particulierement interessant d'elirniner les pixels 

15 correspondant a du bruit en particulier si I'epaisseur moyenne, et done \a% 
zone utile se trouvent a une valeur bien plus elevee (par exemple autouri • 
d'une vingtaine de centimetres) comme ici avec la courbe en tirets T2. ; 

C'est la. raison pour laqueile, sur .la figure 1 le point M a ete regie en * 
tenant compte. de la zone d'anatomie maxiniale. . . - 

20 On voit que pour le cas T2 represente. en tirets sur la figure 4, une!^ 

telle determination aurait conduit a deplacer le point M bien plus vers IgC^ 
droite en M*, c'est a dire a attribuer la dynamique du detecteur a des:M 
epaisseurs bien plus elevees sans perdre une partie importante de cette 
dynamique pour montrer des epaisseurs faibles, dont on aurait determine 

25 par ailleurs qu'elles etaient sans interet puisque representatives uniquement 

de bruit. Sur le plan pratique, d'une maniere presque lineaire le deplacement * 

du point. M vers la droite sera accompagne d'une augmentation de la haute 
tension en kilovolts appliquee entre I'anode et la cathode du tube T. 

Pour mener a bien cette recherche, optionnelle mais preferee, de 

30 I'anatomie maximale, apres I'etape 3 est proposee une etape 4. figure 3. de 
correction de rhistogramme reduit pour tenir compte de la correction 
anatomique max, 

L'etape 4, qui n'est pas indispensable bien que preferee, est suivie 
d'une etape 5 de determination de I'epaisseur moyenne a partir de 
35 rhistogramme reduit, ou de preference a partir de rhistogramme reduit 
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torrige/Gette determination qui sera vue plus loin, comporte un passage du 
domaine des niveaux' de gris au domains des epaisseurs. L'etape 5 est 
suivie par una etape 6 .au cours de laquelle on calcule la dynamique AEPT, 
cette dynamique etant egale, ou correspondant, comme il a ete dit ci-dessus 
5 a recart entre Tepaisseur moyenne EPTmoyen et le seuil d'epaisseur S, ou 
mieux, le seuil d'epaisseur corrige anatomique max. Uetape 6 est suivie 
d'une etape 7 de type connu dans iaquelle on regie Tinstallation en fonction 
. . de I'EPTmoyen ainsi determine et du AERT ainsi calcule. 

Si le praticien ne souhaite pas, pour des raisons qui sont propres a 

10 son metier, beneficier du reglage optimal ainsi obtenu, mais sMl desire 
neanmoins etre debarrasse des zones de saturation, on agit de meme en 
eliminant de Thistogramnie tous les niveaux de gris correspondant 
uniquement a cette saturation. En pratique a titre d'exemple, on elimine tous 
les pixels dont les niveaux de gris correspondent a la saturation, 16384 dans 

15 rexemple, ainsi qu'a ceux qui sont compris dans une bande de niveau de 
gris comportant un certain nombre de niveaux de gris, par exemple 5 
niveaux de gris. Cette bande Bsat est egalement montree sur la figure 1. De 
cette fagon. on ne prendra en compte que les pixels dont les niveaux de gris 
vont de 0 a 16379. On peut alors choisir, tout en gardant la meme epaisseur 

20 moyenne de patient EPTmoyen, puisque celle-ci est exacte, de regler la 
dynamique AEPT en chpisissant un point M" en lieu et place du point M. On 
constate simplement a I'examen de la figure 1, dans ce cas, que si 
I'epaisseur moyenne continue a etre bien calculee, la dynamique n'est pas 
aussi optimale. Elle permet toutefois d'imager les zones B du bord du patient 

25 • sans avoir a supporter d'effet de bord, formant des liseres dans les images, 
et dont la presence peut dans certains cas etre genante. 

En ce qui concerne la determination de I'epaisseur S de seuil en 
dessous de laquelle les contributions des tissus sont jugees sans interet, il 
faut de preference tenir compte de la geometrie d'acquisition de I'installation 

30 au moment de I'acquisition de la premiere image de test a partir de laquelle 
EPTmoyen et AEPT sont determines. A cet effet, I'etape 1 comporte de 
preference une premiere etape 9 au cours de laquelle une valeur AirGap 
d'un espace d'air E est calcule. Cet espace E correspond a I'espace situe 
entre le bord inferieur du corps C du patient et la grille anti-diffusante GA. On 

35 comprend que plus cet espace est grand, plus un phenomene de diffusion 
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Compton. certes marginal mais neanmoins existant, se propagera dans un 
espace plus grand de Timage. Get espace E. de mesure AlrGap, est calcule 
•selon la formule 1 suivante : , 



^ orr. r ^- EPTtheshold 
Air Gap = SID — IsoDis tan ce_ ^ 



Dans cette formule 1 SID represente en cm la distance de la source 
de rayons X, le tube T, a I'jmage, le plan du • detecteur D. La variable 
IsoDistance represente la distance entre le bord inferieur du corps du patient 
et la source de rayonnement X. Dans un exemple correspondant a une 
10 geometrie d'acquisition d'une installation connue. cette variable IsoDistance 
a une valeur de 70,5 cm. En pratique, ces valeurs peuvent etre mesurees 
sur rinstallation utilisee, a moins d'etre disponibles dans des tables dans des 
enregistrements correspondant a des etats d'utilisation de I'lnstaliation-^-^La 
variable EPTthreshold divisee par deux correspond a une valeur purement 
15 arbitraire, typiquement 3 cm, parce que cette valeur est connue dans le 
■ domaine. Une autre valeur que 3 cm peut etre choisie. Elle pourrait 



correspond a Tepaisseur equivalente en dessous de laquelle on est sur 

qu'aucun . tissu interessant n'est present. On a ainsi decouvert que cette 
20 approche de la valeur de Tespace E, independamment de sa significatiori"*en 
^* terme de* longueur, etait particulierement utile pour la prise en compte des 

effets nefastes de la diffusion Compton dans les reglages de rinstallation. 

L'etape 1 comporte ensuite une etape 10 de calcul d'une valeur d'une 

variable ScatterComp, representative de' la diffusion Compton. Cette etape 
25 10 comporte le calcul de la formule 2 suivante : 



.1. 



dependre du lieu d'examen dans le corps C. Cette valeur EPTthreshold "^^^ 



li'. 
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sa 



+ {sbxEPTthrehsold) 
+ (sc X ^SutfaceFdbk x 1 o) \. 
+ {sdx AirGap) 
+ {sexkVp_actual) 
+ {sf X AirGap^) 
+ {sg X EPTthrehsold x ^ SurfaceFdbk x 1 o) 
+ [sh X EPTthrehsold x A irgGap) 
+ {si X EPTthrehsold X kVp _actual) 
+ [sj X ^SwfaceFdbk X 1 0 x AirGap) 
+ [sk X ^SurfaceFdbk x\<^xkVp _actiial] 
+ [si X AirGap x kVp._ actual) . 



Dans cette formule, SurfaceFdbk represente en cm^ la surface du 
champ de vue FOV. et kVp_actual represente la haute tension avec laquelle 
a ete alimente le tube T au moment de I'acquisition de {'image d'essai. Les 
autres variables sont celles vues ci-dessus. Les coefficients sa a si sont des 
coefficients obtenu par une regression. En pratique, on se livre, pour une 
installation donnee a un nombre important de mesures sur des fantdmes de 
densites radiologiques connues, et on determine la variable de diffusion 
Compton ScatterComp par ailleurs. La regression consiste alors a minimiser 
rerreur de determination des douze coefficients sa a si pour le lot 
d'experiences menees. Dans un exemple, ces coefficients ont les valours 
suivantes donnees dans le tableau 1 suivant : 
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oriu Lrase 






sa 


-/ OODD 






SD 


u. 1 oUi^y 1 1 






SC 


-U.U1 UU14^:/i 






sd . 






-U.U/DOOU^O 


se 


u.Uf o/iyooo 






sf 


-2.78306E-05 




0.002540896 


sg 


0.000418987 




0.000775114 


sh 


-0.001951803 




-0.004902506 


si 


-0.002153036 




-0.001676272 


sj 


2.53236E-05 




-2.98979E-05 


sk 


7.9841 3E-05 




6.4485 1E-05 


si 


-0.001169081 




-0.000729517 



Le tableau 1 comporte deux colonnes representant les valeurs des ' 
coefficients sa a si selon qu'une grille GA etait presente, Grid Case, ou non, 
No Grid Case. Les valeurs indiquees dans le tableau 1 ne sent pas uniques > 
5 pour rinvention. Elles dependent entierement de I'installation. Elles peuvent-- 

etre recalculees par regression pour chaque installation. • ' '? 

Rendu a ce stade, il est possible de calculer, selon la formule 3' % 
suivante, correspondant a Tetape 11 de I'etape 1, une variable intermediaire ' . i 
EPTinter: • " ' . " * . ^ 

10 EPTinter = EPTthreshold - ScatterComp \ 

qui est Tepaisseur equivalente, nnaintenant calculee, qui sert de seuil d4 ''i 
partage entre les epaisseurs reputees significatives et cedes qui ne le soinV 
pas. On pourrait se dispenser de ce calcul cbnnplique et imposer une valeur ' 
EPTinter purement arbitraire, qui ne tiendralt en particulier pas compte de la 
15 geometrie de rinstallation et des consequences des diffusions Compton. On 
obtiendrait tlDutefois, sauf hasard improbable, une determination moins 
correcte de I'epaisseur de seuil, et done une determination moins bonne de 
EPTmoyen, et done une sous optimisation du reglage de I'installation. 

La mise en oeuvre de Tetape 2 se realise par le calcul de la formule 4 
20 suivante : 



SFBthreshold = exp(^?i2 -{-Wii - i - /« + bii)) 



dans laquelle <1> est une fonction dite de Tansig, connue, et est donnee par 
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• la formule 5 suivante .: , ' . " ' 

Dans la formule 4, le resultat SFBthreshold correspond a I'equivalent 
en dose de I'epaisseur seuil S qu'on a voulu retenir pour eiiminer les zones 
5 de. non-anatomie. On effectue .au bespin par ailleurs une autre 
. transformation de SFBthreshold. de type connu. pour passer du seuil en 
. dose a un seuil en niveau de gris. Mais il est toutefois possible de travailler 
dirictement en dose, sans passer par les niveaux de gris. Les termes b de la 
formuie 4 sont des vecteurs, et les termes W sont de matrices. Les. 
10 dimensions de ces vecteurs et de ces matrices sont donn§es par le tableau 
' :2 suivant, qui comporte, dans les memos conditions que le tableau 1, 
precedent uh cas ou une grille antidiffusante est utillsee, Grid Case, et un 
cas ou une telle grille ne i'est pas, No Grid case. Les valeurs indiquees dans 
le tableau 2 ne sont pas uniques pour I'invention. Elles dependent 
15 entierement de I'installation. Elles peuvent etre recalculees par regression 
pour chaque installation . • 
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column4 
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-1.352894e+000 


column4 


0.2825646 




1.5862096+001 


column4 


-0.2311099 




8.932972e-002 


column4 


-0.05804851 




-2.6045966-001 


columns 


- 0.004109868 




-1.3884686-002 


columns 


0.00886492 




1.4581866+001 


columns 


-0.007394629 




3.8668626-003 


columns 


-0.001900335 




-9.5012906-003 



Le terme In de la formule 4 est un vecteur a 5 dimensions donne par 



rexpression 1 suivahte: 



In = [Patient _ size _ normal, kVpnonnal^ mAnormal, cu _ thickness normal, al _ thichtess _ normal]^ 



5 Dans cette expression 1 Patient_size_normal= 0.6 * EPTinter. En outre 
chaque valeur "normal" est de preference une valeur normalisee 
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correspondant :aux formules suivantes : 
kVp normal est donne par 



kVpnormal ■ 



kVp _ actual -^kVp min 



kVp max— kVp min 



Dans cette. formule kVpmin et Wpmax sont Les valeurs. minimales et 
maximales de la tension d'utilisation, alors que kVp_actual represente la 
haute tension de I'installatibn dans les conditions de I'image test. 
Dans cette expression 1 mAnornnal est donne par: 



mA normal = 



mAs actual x 



' SID ' . 



mr 



mas 



mr ___mas _cal j 



Dans cette formule SID est la valeur de Source - to - Distance, SID_EPTnn 
10 vaut 100 cm. mR_mAs_cal = 4.0858, et mR^mAs est la valeur catibree de 

mR/mAs " 

Dans cette expresssion 1 Cu_thickness_normal et 
ALthickness_normal sont donnes par la formule suivante: 



filternoifnalQi) = 



Spectral„Filter_feedback(thickness[ 1 ] - Spectral_filter__min[l ] 
SpecrtalFilterMaxThichness[\] - Spectral_filter__min[l] 



15 Dans laquelle les valeurs des filtres Cu_lo, Cu__hi, ALIo et ALhi sont 

respectivement les epaisseurs minimales et maximales de filtration en cuivre 
et en aluminium de Tinstallation, En appliquant la formule 4, on montre done 
qu'on peut realiser I'etape 2 c'est-a-dire calculer uniquement a partir des 
parametres de reglage de Tinstallation, le niveau de gris correspondant au 

20 seuil 8. En pratique, ce calcul peut etre mene avant meme que I'image ne 
soit acquise et que son traitement selon les etapes 3 et suivantes ne soit 
entrepris pour editer les valeurs EPTmoyen et AEPT recherchees. 

La question est aussi de savoir comment transposer Tenseignement 
ci-dessus selon I'invention a une autre installation. Cette transposition est 

25 rendue possible en modelisant la formule 2 pour une installation donnee, 
Une telle modelisation peut etre menee a Taide d'un reseau de neurones. 
Plus theoriquement, une telle modelisation est realisee en menant un grand 
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nombre d'experiences. Pour ces experiences, on fait varier les parametres 
d'entree In et on mesure les doses regues en correspondance. Puis on 
minimise I'erreur de calcul de la correspondance entre ces doses et ces 
parametres d'entree, selon cette formule 2, en recherchant les composantes 
5 des termes b et W qui satisfont le mieux cette minimisation. Le modele de la 
formule 4 est connu par ailleurs, La particularite de Tinvention revient a 
Tutiliser en inverse, • 

De. preference, la determination des parametres optinnaux de 
Tacquisition de radiographie, ainsi que le reglage de Tinstallatioa, sont 
10 realises en temps reel. Une fois acquis prealablement les coefficients des 
tableaux 1 et 2, pour une installation donnee, les calculs des conditions 
optimales de reglage de Pinstallation sont realises en quelques millisecondes 
apres acquisition de Timage test. 

"r 
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REVENDICATIONS 

1 - Procede de determination des parametres "optimaux d'une 
aequisition de radiographie, dans lequel, 
5 - on acquiert une premiere image test du patient dans des conditions 

de reglage connues d'une installation de radiologie utilisee a cat effet, 

- on mesure pour ces conditions de reglage connues une epaisseur 
**moyenne du patient a partir de cette premiere image test, et 

V. ' - on determine les parametres optimaux de Tacquisition a partir de 
10 cette epaisseur moyenne, 

caracterise en ce que 

- pour mesurer ('epaisseur moyenne a partir de cette premiere image 
test, on ecarte dans cette premiere image test des elements d'images qui ne 

. ^ sont pas representatifs de parties anatomiques significatives du patient. 
15 2 - Procede selon la revendication l, caracterise en ce que 

- on regie Tinstallation de radiologie en fonction de ces parametres 
optimaux, 

- notamment on regie la haute tension appliquee entre une anode et 
une cathode d'un tube a rayon X de cette installation, ce reglage etant 

20 fonction de I'epaisseur moyenne du patient examine 

. 3 - Procede selon Tune des revendications 1 a 2, caracterise en ce 

que 

- pour ecarter, on elimine de ia premiere image des elements d'image 
dont une caracteristique dans I'image est situee au-dela d'un seuil (S), 

25 - ce seuil correspondant a une epaisseur de limite d'interet de corps 

de patient, 

- en dessous de cette epaisseur de limite d'interet, on considers que 
rimage n'a pas d'interet, 

- et on calcule I'epaisseur moyenne a partir d'uri histogramme reduit 
30 d'elements d'image, dans cet histogramme reduit les elements d'images 

elimines ne sont pas presents. 

4 - Procede selon la revendication 3, caracterise en ce que 

- on prend comme valeur d'epaisseur moyenne la moyenne des 
epaisseurs equivalentes correspondant aux elements d'images de la 

35 population d'elements d'image de Thistogramnie reduit. 
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5 - Procede selon la revendication 4, caracterise en ce que 

- on soustrait de la population d'elements dimage de rhistogramme 
reduit un nombre donne d'elements d*image, 

- les elements damages soustraits sont ceux dont les epaisseurs 
equivalentes sont les plus faibles a partir du seuil/ 

6 - Procede selon Tune des rievendications 3 a 5, caracterise en ce 
que pour le reglage de rinstallation, 

- on retrouve par analyse inverse, a partir des parametres de reglage 
de rinstallation pour {'image test, un seuil d'element d'image, par exemple un 
seuil en niveaux de gris ou en doses, qui correspond a une epaisseur au- 
dela de laquelle les zones de tissus sont reputees ne pas avoir d'interet. 

7 - Procede selon la revendication 6, caracterise en ce que, avant 
I'analyse inverse, 

- on corrige I'epaisseur en^ fonction d'une epaisseur arbitraire, d'une 
geometrie d'acquisition de Timage et en fonction d'un phenomene dej 
diffusion Compton qui en resulte selon les formules y. 



A irGap = SID - IsoDis tan ce 



EPTthreshold 
2 



et 



sa 



+ [sb X EPTthrehsold) 




ScatterComp 



^{sfxAirGap^) 




8 - Procede selon Tune des revendications 6 a 7, caracterise en ce 
que, pour Tanaiyse inverse, 



1 er depot'^ 



22 : 

- on calcule le resultat de la formule 

SFBthreshold = exp(b\2 + Wn- ^(W\ \'In + b\ i)) 

dans laquelle, In represente les conditions de reglage connues de 
rinstallation, les bij et Wij sont respectivement des vecteurs et des matrices, 

- on calcule les bij et Wij par apprentissage, notamment par. 
5 minimisation d'une erreur de caicul de cette formule pour un ensemble de 

seuils, SFB threshold, et pour un ensemble de conditions variees de reglage 
derinstallation. 

9 - Precede selon Tune des revendications 3 a 8, caracterise en ce 
que pour regler rinstallation, 

10 - on cale la dynamique d'un detecteur de ('installation de telle fagon 

qu'une proportion donnee d'un maximum de cette dynamique corresponde a 
une tranche d'epaisseur equivalente, cette tranche etant comprise entre 
Tepaisseur moyenne du patient et une epaisseur correspondant au seuil. 

10 - Precede selon Tune des revendications 3 a 8, caracterise en ce 
15 que pour regler installation, 

- on cale la dynamique d'un detecteur de rinstallation de telle fagon 
que cette dynamique corresponde a une tranche d'epaisseur equivalente, 
cette tranche etant comprise entre Tepaisseur moyenne du patient et une 
epaisseur du patient la plus fine encore visible dans I'image avant saturation. 
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Anatomic 
max 



noir 0 



Nombre de 
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Calcul de 
Air gap 



Calcul de 
ScQttercomp 



Calcul de 
EPT inter 



.10 
,1 
11 



Calcul par analyse inverse 
de SFBth seuil et gth " 



Realisation de rhistogramme 
et de I ' histogramme reduit " 



Correction de 
I 'histogramme reduit"^ 



A 



Determination de EPT moyen 
a partir de rhistogramme ^ 
reduit corrige 



Calcul de la 
dynamique AEPT' 



Reglage en fonction 
deEPTm et A EPT" 



Reglage en fonction , 
de EPT et global max 



,8 
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Zsat 



Nombre de Epa'^^f Epaisseur 
pixels A mimmale 2 mmimale 1 



Global max 




Niveoux de gris 
du seuil d' epaisseur 



Niveaux 
de gris 



Fog. 4 
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